INP MSU Preprint-95-2/366 
Journal-ref: Int. J. of Mod. Phys. C7, No.6 (1996) pp.761-774 

18 December 1996 



RECURSIVE ALGORITHM FOR THE GENERATION OF RELATIVISTIC 
KINEMATICS FOR COLLISIONS AND DECAYS WITH REGULARIZATIONS 

OF SHARP PEAKS 



V.A. Ilyin, D.N. Kovalenko, A.E. Pukhov 

Institute of Nuclear Physics of Moscow State University, 119899 Moscow, Russia 



Abstract 

We describe an algorithm for the generation of relativistic kinematics for collision and 
decay processes with multiparticle final states. In the framework of this algorithm 
it is possible to generate different kinematics covering most of practically interesting 
cases. One gets a possibility to introduce different sets of integration variables. As a 
result different sets of kinematical singularities can be regularized. To smooth sharp 
peaks some regularization formulas and procedures are used covering most typical 
cases. The algorithm is realized in the package CompHEP created for automatic 
calculations of collision and decay processes. 



1 Introduction 



To evaluate numerically cross sections and other characteristics of elementary particle collisions and to calculate 
decays of unstable particles it is necessary to integrate squared matrix element over the phase space. 

The first step is to express the integrand in terms of independent kinematical variables — integration 
variables. Note that after evaluation of Feynman diagrams the integrand is expressed in terms of scalar products 
of Lorentz momenta of in- and out- particles, and polarization vectors. 

The next step is connected with a singular behavior of the integrand that is typical case in modern col- 
lider physics calculations (see Ref. [Q, for example). In many cases straightforward methods of the numerical 
integration fail because of sharp peaks. Standard procedure to treat carefully such peaks is to transform in- 
tegration variables in order to get nonsingular function in the new variables. We will call such procedure the 
regularization. 

These two steps are closely connected. To make it clear let us note that singularities appear as zeros (more 
exactly — almost zeros) of propagators denominators. So one can say that each singularity is connected di- 
rectly with a squared Lorentz 4- vector expressed in its order in terms of linear combination of particle momenta. 
However it is clear that singularity can be smoothed by integration variable transformation only if the corre- 
sponding invariant is one of these variables or related directly to such variable. Thus one can integrate carefully 
if both steps will be performed in proper relation to the set of singularities. This problem becomes nontrivial 
especially in so-called automatic calculation of collision and decay processes. This new approach, based on 
high perfomance and user-friendly programs, is developed now intensively Refs. 0-0- With the help of these 
programs one can calculate different processes in the framework of the various physical models. It is clear that 
the generation of kinematics is an important part of this approach. 

There are many known kinematical schemes implemented in existing approaches (see Ref. [gj , for example) . 
In GRACE || special library of kinematics, corresponded to different sets of kinematical singularities, was 
created. In EXCALIBUR || a multichannel appproach to compute four fermion processes in e + e~ annihilation 
is presented. In this paper we demonstrate that the problem discussed can be treated in the framework of some 
general algorithm based on the recursive representation of kinematics as a tree of elementary decays 1 — > 2. It 
opens the possibility to generate kinematics (choose integration variables, express squared matrix elements in 
terms of these variables and introduce proper regularizations) for almost all kinds of processes with the help of 
one program with a high level of user's interface. This algorithm was realized in the package CompHEP |2j and 
representative computational examples are given here. 



2 Kinematics 

By definition, the kinematical scheme is a way of phase space parameterization. 

We define the kinematical scheme as a tree of elementary decays. Under the elementary decay we understand 
the parameterization in 2-particle phase space. Note that the words elementary decay do not correspond to real 
1^2 decay in most cases. We use this terminology for convenience. 

Below we consider multiparticle phase space for collisions. For decays the corresponding formulas can be 
written in a similar way. 

Let us take any two out-particles with the numbers i and j. We consider their momenta as momenta 
produced in the following elementary decay: 

Qij -> Qi , Qj , (1) 
where qij = qi + qj = (Eij, qij). Let us introduce a parameterization for this elementary decay in the rest frame 
of decayed 4- vector as 

R ' 

dT ^l^ (2^)3.4^ -^ - dQi - (2) 

Here the measure dfli corresponds to the angular part of dqi. Then Ri = \qi\ in the rest frame of q^, and 
Sij = qfj- 

We repeat now this construction recursively: on each step two cluster 4- vectors qi,qj are taken from the 
rest part of out-momenta or from the decayed 4- vectors qu appeared on previous steps. After (N — 1) steps we 
get the following formula for differential cross section (|M| 2 is a squared matrix element — the normalization 
corresponds to the conventions 0): 

N-l 

do = 7p^\M\ 2 ■ dT* 2 {h, Ji) • [] dT 2 (Ii, J{) , (3) 
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where the 2-particle phase space measure for the 1st elementary decay p\ + P2 ^ <hi i <7,/i does not include 
integration over s = {pi + P2) 2 (here pi and P2 are momenta of colliding particles): 

dTUiuJ^j^^-dn,. (4) 

Here Ri 1 = {IgjJ, in CMS}, by CMS we denote the kinematical frame of center mass system. 

The angular measure dfli can be written in terms of two angular coordinates of the 3- vector <f/. Let us 
introduce two reference 4- vectors Vi — (Vj , Pi) and Ai — {A®, Ai) for some elementary decay in the rest frame 
of the decayed vector qjj. So, let the polar coordinate of <fj, the angle 9/, be counted from the 3-vector Vj. 
Then let us take the azimuth coordinate of qj as the angle between two planes (Pi, Ai) and (Pi, qi). As a 
result the measure dfij takes the form 

dfij = dcos6/ • d$j . (5) 

Taking the reference vectors as different combination of the momenta of in- and out-particles we can connect 
different Lorentz invariants with introduced integration variables. Let us consider, for example, elementary 
decay qu —* qj, qj and Vi — J2 n Pni where p n are momenta of some in- and out-particles. The relation between 
the angle 0/ and the invariant {^2p n + Qi) 2 can be written in the following form: 

(]>>« + ?/) 2 = (J2p^T + ^ 2 + 2 {T,Pn) ■ i°i- 2 |EH-i^i cose '- ( 6 ) 

In the following we use this relation to introduce regularizations of sharp peaksQ 
Let us introduce numbering of the elementary decays by the following rules: 

1 . initial elementary decay (p\ + P2 — * qi t , qj t ) has a number i = 1 ; 

2. all elementary decays on branches, growing from the clusters U and Ji of the i-th elementary decay, have 
numbers larger than i; 

3. any elementary decay on the branch growing from the cluster Jj has a number smaller than any elementary 
decay on the branch growing from the second cluster Ji of the i-th elementary decay. 

We will use a terminology 'previous elementary decay'' if it has smaller number than the current one. In Com- 
pHEP realization the integrations in (3) are ordered in the correspondence with this numbering of elementary 
decays. 

2.1 Kinematics summary 

• The kinematical scheme is a tree of "elementary decays". 

• The number of elementary decays equals (A^ out — 1), where JV out is a number of out-particles. 

• Each elementary decay is defined by 

— three 4-vectors: decayed qu and two cluster vectors, qi and qj, with qjj = qj + qj; 

— reference vector Vi which can be chosen from some combinations of particle momenta. These com- 
binations have to be already defined as cluster momenta or decayed vectors in previous elementary 
decays. 

• The vectors qu and Vi has to be linearly independent. 

• Each elementary decay qu — * qi, qj adds three new integration variables: 

— su — squared decayed vector qu', 

— Qi and — polar and azimuth angles of qi counted from the reference 3-vectors Vj and Ai in the 
rest frame of qu- For decays or collisions with fixed energies of in-particles the 1st elementary decay 
gives only two independent variables, 0i and $1. In that case the number of integration variables 
equals 3 • (N out - 1) - 1 = 3 • iV out - 4. 

x The relation between angle $j and some Lorentz invariant is more complicated because it includes also cos0. Of course, it 
is possible to use the angle for the regularization purposes. However, we tested such regularizations and found them to be not 
effective and we do not consider this possibility here. For integration over $7 we take some fixed reference vector Ai- It could be 
unit vector orthogonal to the collision axes, for example. 
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3 Invariant Cuts 

In the framework of the algorithm described above it is possible to introduce effectively kinematical cuts over 
some combinations of Lorentz momenta of in- and out-particlcs. 

In general we can introduce an invariant cut in two ways. The first is more effective and it can be realized 
if the corresponding invariant is connected directly with one of the integration variables of kinematical scheme 
(e.g., as Eq. ^)). In this case the corresponding cut can be implemented by the limitation of the integration 
volume. 

In other cases entered cuts can be applied by multiplication of squared matrix element on the corresponding 
step function. This way can be not so effective, especially in the case of strong cut when the integrand equals 
zero in almost all phase space volume. 



4 Regular izat ion 

In some cases the squared matrix element has singularities when the integrand has sharp peak(s) in phase space 
and the standard Monte Carlo integrator (even adaptive one) could fail. The typical examples are given below. 

4.1 General regularization formula 

Let the integrand F(x) to be represented as F(x) = f(x) ■ (gi(x) + • ■ ■ + g n (x)), where f(x) is a smooth function 
and functions gi{x) have sharp peaks. One can transform the integration variable in the following way: 

1 " 
dx = — dy, y{x) =^2,Gi{x) = G{x) , (7) 

X>(*(V)) 
i=l 

where integral functions are Gi(x) — f x gi(x')dx' . The function x(y) = G~ 1 (y) we will find from the relation 
J27=i Gi{x) — y = numerically. Also let us transform the integration volume to the interval [0, 1]: dy — 
(G(b) - G(a))dy, 0<y<l. 

As a result the integral can be rewritten as 



^9i{x{v)) )J{g}(y)dy, (8) 




where Jacobian equals J{ g }{y) = — g.( g ( y )) (^(ft) — ^( a ))- 

We see that the singular integral is transformed now to smooth integration of the function f(x). Of course 
integrand can have more complicated functional structure then we have considered. However, if we know the 
structure of singularities it is not so difficult to write some model function of the form (gi (x) + ■ ■ -+g n (x)) which 
will match each singularity. And as a result the transformation of the integration variable described above will 
regularize the integrand. Here we would like to stress again that the structure of singularities can be derived 
easily from the analysis of denominators of calculated Feynman diagrams. 

Note that it is possible to smooth several peaks over one variable appearing from the propagators of different 
Feynman diagrams. In this case the model singular function gt(x) is constructed with gi(x) corresponding 

to separate peaks. 

Then we consider the constant g n (x) — 1 as one of model function. It is necessary to avoid neglect of 
contributions of nonsingular areas after regularization. 

4.2 Types of regularizations 
4.2.1 ^-channel pole 

Let us consider a typical example when in some Feynman diagram the in-particle with momentum pi radiates the 
virtual particle with mass m and after that transforms to the out-particle with momentum p 3 . The corresponding 
propagator of the virtual particle is 

r^-j, t^( P3 ~ Pl ) 2 . (9) 
t — m z 
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Figure 1: Feynman diagram of7 + e — > e + Z + H with i-channel type of singularity, produced by the 
exchange of virtual electron. 
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Figure 2: Feynman diagrams of e + 7 — * e + W + + W with ^-channel type of singularity, produced by the 
exchange of virtual photon. 

The integrand will have sharp peak in t near its upper bound t max if m 2 and i ma x are much less than the energy 
of in-particles: 

m 2 < s , i max < s . (10) 
Typical example is given in Fig. |l| by the Feynman diagram of the process 

7(pi) + e~(p 2 ) -> e-(ps) + Z(p 4 ) + H(p 8 ) , 

where m = m e is the electron mass, and we indicated momenta of particles in the brackets. 
Another example is given by the process 

e~(pi) + 7(R2) -» e_ (P3) + ^+(p 4 ) + W~(P6) , 

when m-electron radiates the virtual photon (here m = 0). The Feynman diagrams for this case are shown in 
Fig-H 

Separate squared diagrams can include terms asymptotically proportional to 1/t 2 and also to 1/t. The 1/t 2 
terms coming from different diagrams can cancel each other up to the terms of 1/t type (in the m e — ► limit) 
due to the gauge invariance. However, for processes 2 — > n, with n > 2, if the singular exchange propagator has 
its momentum composed from more than one out-particle momenta then there is no cancellation of 1/t 2 terms. 
For some processes of 2 — ► 2 type, like elastic scattering of two charged particles, there are no cancellation of 
1/t 2 terms also. In these cases we have much stronger singularity. 

Let us consider a i-channel singularity over tj = (p m — Y^, Pi Ut ) 2 i where p ln is momentum of one of in-particles 
and sum is made over the some subset of out-particles. Suppose that we constructed a tree of elementary decays 
in such a way that we can choose some elementary decay, with clusters momenta q\ and qi , and the reference 
vector V, as a combination of clusters momenta from previous elementary decays, according to one of two 
variants: 

(p in -$>r) 2 '' = (qi+V) 2 , or [r-Y^PT') 2 = ^-2+V) 2 . (11) 
In these cases we have 

tj = (q) 2 + (V) 2 + 2(q V -\q\\P\ cos 6) , q = q 1>2 . (12) 

So, the invariant tj will be linear function of the integration variable cosO. 
In the case of cancellation of 1 /t 2 terms the model singular function is 

^° S9 )^- (13) 



4 




Figure 3: Feynman diagram with decayed W propagator contributed in the process e + e + — > + /i + W + . 
Its integral function is 

Gi (cos G) = In | cos6 -C |, (14) 

where 

= m 2 - (V) 2 - (q) 2 + 2qoV 
2\V\\q\ 

Here m 2 is the squared mass of the corresponding singular propagator (Eq. (^|)). 
In the cases without cancellation of 1 /t 2 terms the model singular function is 

gi(cos9) = (cose- c r (15) 

with integral function 

Gi(cos6) = - 1 (16) 
Go — cos H 

4.2.2 s-channel decay pole 

When a process includes the s-channel decay of some unstable particle it is necessary to introduce a finite width 
in the corresponding propagator to get a physically meaningful (and finite) cross section. Such diagrams will 
give in squared matrix element the factor 

(s fc - M 2 ) 2 + T 2 ■ M 2 ' 

where M is the mass of unstable particle, T is its width and Sk is the squared invariant mass of decay products. 
The integrand has a sharp peak if 

r • M < a . (18) 

An example is the process 

e" + e+ -> Df, + [i + W + 

which includes the W~ decay into the lepton pair v^^i. The Feynman diagram with such singularity is shown 
in Fig. |. 

If Sk is an integration variable (so it is squared decayed momentum of one elementary decay of the kinematical 
scheme) the corresponding model singular function can be taken as 

9M = {3k _M, )2 + £2 (19) 

and its integral function equals 

Gi(s k ) = X - arctan f Sfc ~ M \ , (20) 

where e = T ■ M. 

Let us now consider the case when in some kinematical scheme Sk does not correspond to decayed momentum 
of any elementary decay. So we can not regularize this singularity by the previous way. To carry out the 

regularization in this case we need to choose elementary decay, with that 4-momenta qi and q 2 of its two 
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Figure 4: Feynman diagrams for e +e + ^e +e + + Z with infrared type of kinematical singularity. 

clusters, and the reference vector V , as a combination of clusters momenta from previous elementary decays, 
according to one of the following variants: 

^■=(E^ Ut ) 2 = (<Zi+7>) 2 , or s k =(j2P° Ut ) 2 = (Q2+V) 2 . (21) 
In these cases we have 

s k = (q) 2 + (V) 2 + 2(q V -\q\ \V\ cos 6) , (7 = 91,2- (22) 

Thus the invariant Sk will be linear function of the integration variable cos0. The model singular function for 
such cases is 

&(cos6) = 7 — , 2 (23) 

(cos 6 - C/ ) + e o 



and its integral function equals 



where 



and 



M 2 - (q) 2 - (V) 2 - 2q V 
-2\q\\P\ 

e 



Gi ( cos 9) = — arctan f C ° s9 C ° I , (24) 



C = 



4.2.3 Infrared type of the singularity 

Let two light out-particles with masses mj and m 2 (for instance photon, electron, u or d quarks) be "products 
of the decay" of some virtual particle with mass m, and m < m\ + ni2- Sharp peak arises from the propagator 
of the "decaying" virtual particle with invariant mass Sf. 

' Si > (mi + m 2 ) 2 , (25) 



(sj — m 2 ) 2 
if the masses mi, Tri2, and m are small: 

mi, 77i2, m <C • (26) 

Note that in the Standard Model such singular structure appears typically if one of these three particles is 
photon or gluon. An example is a "decay" of virtual photon into two light charged particles. Such peak arises 
in some Feynman diagrams (see Fig. @) of the process 



In this case m = and mi = 7772 = 77i e — mass of electron. 

The singularity of infrared type is present also in the case of out-photon emission from the charged out- 
particle. In this case mi = and m 2 = 777 is a mass of charged particle. The total cross section is infinite in 
this case — there is a nonintegrable pole at zero photon energy (so-called infrared catastrophe). However one 
can introduce the cut 

Si > S 2 +m 2 , St = (p mi +p m2 ) 2 , (27) 
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where p mi and p m2 are momenta of the out-photon and charged out-particle. In this case sharp peak appears if 

5<V^- (28) 

All these situations are different from the previous example of s-channel singularity. First, the pole is placed 
beyond the physical region. Second, due to the gauge invariance, the singularity will be of the first order, terms 
of l/(sj — to 2 ) 2 type will always cancel each other up to the terms of l/(s<i — to 2 ) type. So we have the similar 
functional problem as in the case of i-channel pole and we can use the formulas of the f-channel regularization 
to handle such poles. 

5 Computational Examples 

The method described in the previous sections was realized in CompHEP Q package. Let us inspect the 
efficiency of this method in some examples. The results of cross section calculations by CompHEP with Monte 
Carlo integration are listed below. The program VEGAS S was used for multidimensional adaptive MC 
integration. 

5.1 i-channel singularity 

Let us begin from the process 

e_ (Pi) +l(P2) -> e~(p 3 ) + Z(p 4 ) +H(p 5 ) 

with singularity appearing from the exchange of electron in f-channel (Fig. ^). First let us use the default 
kinematics in CompHEP: 

1st elementary decay: p\ + P2 — * P3, Pa + P5 with V = p\\ 
2nd elementary decay: p 4 + p§ — ► P4, ps with V = —p\ + p$ . 

For the number of points per iteration equals 10000, we can see the following CompHEP screen output for the 
total cross section calculation at s/s = 500 GeV: 



IT 


CURRENT RESULTS 
Cross section [pb] 


Error '/, 


ACCUMULATED 
Cross section 


RESULTS 
[pb] Error '/ 


CHI**2 


1 


3. 


685360E-04 


24. 


.68 


3 


.685360E-04 


24. 


.68 


.00 


2 


6. 


. 255344E-04 


27. 


. 15 


4 


.258167E-04 


18. 


.83 


1.78 


3 


4. 


.638301E-04 


10. 


.09 


4 


.541736E-04 


8. 


.90 


.97 


4 


5. 


. 303054E-04 


6 


.19 


5 


.000283E-04 


5. 


.10 


1.36 


5 


6. 


214947E-04 


16. 


.34 


5 


.072243E-04 


4. 


.87 


1.36 


6 


6. 


. 355495E-04 


8 


.07 


5 


.314256E-04 


4. 


.19 


2.10 


7 


7. 


. 45841 1E-04 


7. 


. 14 


5 


.633162E-04 


3. 


.65 


4.05 


8 


9. 


. 273354E-04 


19. 


.66 


5 


.678777E-04 


3. 


.59 


4.04 


9 


8. 


410235E-04 


6 


,03 


6 


.059285E-04 


3. 


.13 


6.65 


10 


7. 


. 566958E-04 


1. 


.20 


7 


.285007E-04 


1. 


.12 


11.64 



Here IT is iteration number, CURRENT RESULTS are results of the integration for current iteration. 
ACCUMULATED RESULTS are cumulative value of the integration of current and previous iterations. These 
numbers , as well as the error and parameter x 2 , are evaluated according standard formulas (see details in 
Ref. §). 

The result above shows that there is no stability from iteration to iteration (see 2nd column — values of 
current results; and 3rd column — current errors). Also there is no reliability with due to large value of \ 2 
parameter. 

Let us now introduce the kincmatical regularization of the t-channel type for the invariant (p^ — P2) 2 with 
mass to = rn e . Used kincmatical scheme allows to introduce this regularization because V in the 1st elementary 
decay equals p\ and p\ = —pi in the 1st decay. The CompHEP screen output of the integration with entered 
regularization will be: 



IT 


CURRENT RESULTS 


ACCUMULATED RESULTS 


CHI**2 




Cross section [pb] Error '/, 


Cross section [pb] Error 7, 
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1 7.621777E-04 .33 7.621777E-04 .33 .00 

2 7.598378E-04 .21 7.605149E-04 .18 .62 

3 7.590317E-04 .18 7.597725E-04 .13 .61 



5.2 s-channel singularity 

Let us consider an example with s-channel singularity (see Fig. ||): 

e~(pi) + e+(p 2 ) i/„(p 3 ) + /i(p 4 ) + ^ + 0p 5 ) • 

If we use the following kinematical scheme 

1st elementary decay: p\ + pi ™ > P3, Pi + P5 with V = pi ; 
2nd elementary decay: pi + P5 ~* Pi, Pb with V = P3 , 

then the CompHEP screen output with 25000 points per iteration and s/s = 500 GeV is: 



IT 


CURRENT RESULTS 
Cross section [pb] 


Error '/, 


ACCUMULATED 
Cross section 


RESULTS 
[pb] Error V, 


CHI**2 


1 


6.407551E-01 


27. 


.60 


6 


.407551E-01 


27. 


.60 


.00 


2 


8.139560E-01 


10. 


.75 


7 


.798860E-01 


10. 


.06 


.77 


3 


7.667456E-01 


3 


.48 


7 


.681075E-01 


3. 


.29 


.40 


4 


8.065722E-01 


2 


,70 


7 


.901984E-01 


2. 


.08 


.71 


5 


8.483987E-01 


4 


.06 


8 


.010476E-01 


1. 


.86 


1.11 


6 


8.442694E-01 


2 


.95 


8 


. 123698E-01 


1. 


.57 


1.33 


7 


8.276117E-01 


2 


.50 


8 


. 165697E-01 


1. 


.33 


1.18 


8 


8.402845E-01 


2 


,54 


8 


.214562E-01 


1. 


.18 


1.15 


9 


8.512594E-01 


3 


.52 


8 


.242752E-01 


1. 


.12 


1.12 


10 


8.232579E-01 


3 


.04 


8 


.241540E-01 


1. 


.05 


.99 



We see that the error of current iteration (see 3rd column) is standing at the level about 3% after 10 iterations 
that it is not satisfactory in many applications. 
Let us choose other kinematical scheme 

1st elementary decay: p\ + P2 — ► P3 + Pi, Pb with V =Pa 
2nd elementary decay: p 3 + p± — > p$, p^ with V = p\ . 



One can see much better results even without regularization: 



IT 


CURRENT RESULTS 
Cross section [pb] 


Error 


ACCUMULATED 
Cross section 


RESULTS 
[pb] Error '/, 


CHI**2 


1 


9.342104E-01 


29.08 


9.342104E-01 


29.08 


.00 


2 


8.762828E-01 


8.78 


8.805859E-01 


8.41 


.04 


3 


8.677625E-01 


1.18 


8.680033E-01 


1.17 


.04 


4 


8.463662E-01 


.65 


8.513031E-01 


.57 


1.19 


5 


8.471336E-01 


.57 


8.492184E-01 


.40 


.99 



Then introducing regularization over invariant mass (p3 +P4) 2 with the mass parameter M = Myy (T^-boson 
mass) and the width parameter T = T\y (M^-boson width) we get more effective integration than in the previous 
case: 



IT 


CURRENT RESULTS 




ACCUMULATED RESULTS 


CHI**2 




Cross section [pb] 


Error '/, 


Cross section [pb] Error '/, 




1 


8.535685E-01 


4.72 


8.535685E-01 4.72 


.00 


2 


8.539195E-01 


.94 


8.539060E-01 .93 


.00 


3 


8.544500E-01 


.60 


8.542897E-01 .50 


.00 


4 


8.449297E-01 


.58 


8.502279E-01 .38 


.69 


5 


8.481124E-01 


.57 


8.495686E-01 .32 


.55 
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The result of regularization is strong decrease of the error of first iterations. 

However, it can happen that for some reason one cannot choose kinematical scheme where this s-channel 
singularity is related to some integration variable as an invariant mass of elementary decay. In the process 
e~ + e + — * e~ + 9 e + W + , for example, when discussed s-channel singularity contributes, other singularities 
are present and they can force us to introduce another kinematics to regularize ihese peaks. In such cases 
one can use regularization of the s-channel singularity using polar angle (Eq. fl22|)). In our demonstration 
calculation let us return to the first kinematical scheme given at the beginning of this section and introduce the 
regularization for the corresponding (p^ + P4) 2 invariant. It is possible to do so because p% has been choosen in 
this kinematical scheme as reference vector P in the 2nd elementary decay. The CompHEP screen output for 
this case will be: 



IT 


CURRENT RESULTS 
Cross section [pb] 


Error '/, 


ACCUMULATED 
Cross section 


RESULTS 
[pb] Error V, 


CHI**2 


1 


8.854719E-01 


4.85 


8.854719E-01 


4.85 


.00 


2 


8.453821E-01 


1.00 


8.468828E-01 


.98 


.84 


3 


8.477949E-01 


.77 


8.474495E-01 


.60 


.42 


4 


8.497170E-01 


.82 


8.482457E-01 


.49 


.30 


5 


8.503980E-01 


.75 


8.488813E-01 


.41 


.25 



One has to compare this screen output with the first one at the beginning of this section. We see that 
the integration is stable with fast convergency, smaller error and reliable value of x 2 parameter, similar to the 
previous case when the regularization over (p^ + P4) 2 as invariant mass of elementary decay was performed. 



Acknowledgments 

This work was partially supported by European Association INTAS (grant 93-1180-ext), Russian Foundation 
for Basic Reseach (RFBR grant 96-02-19773a) and JSPS (Japan Society for Promotion Science). 

We are grateful to E. E. Boos, M. N. Dubinin, and S. A. Shichanin for fruitful discussions and testing the 
program, and to Y. Shimizu for useful discussions. 



References 

[1] D. Bardin et ai, "WW physics and generators," in Physics at LEP II, ed. G. Altarelli, T. Sjostrand, and 
F. Zwirner, CERN report 96-01, 1996. 

[2] E. Boos et ai, SNUTP 94-116 Seoul 1994, INP MSU-94-36/358 Moscow 1994; jhcp-ph/9503280| . 

[3] GRACE manual, KEK Report 92-19, 1993. 

[4] J. Kiiblbeck, H. Eck, and R. Mertig, "Computeralgebraic generation and calculation of Feynman graphs 
using FeynArts and FeynCalc," in Proc. of the Second International Workshop on Software Engineering, 
Artificial Intelligence and Expert Systems in High Energy and Nuclear Physics, ed. D. Perret-Gallix (1992). 

[5] E. Byckling and K. Kajantie, Particle Kinematics (John Wiley & Sons, London, 1973). 

[6] R. Kleiss and R. Pittau, "Weight optimization in multichannel Monte Carlo," NIKHEF-H 94-17 and 
INLO-PUB-4/94; F.A. Berends, R. Kleiss, and R. Pittau, "EXCALIBUR — a Monte Carlo program to 
evaluate all four fermion processes at LEP 200 and beyond," INLO-PUB-12/94. 

[7] "Particle Data Group," Phys. Rev. D50, No. 3 (1994). 

[8] G.P. Lepage, J. Comp. Phys. 27, 192 (1978). 



9 



